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^ ! Abstract 

fSJ ! Relativistic mean field theory with the NL3 force is used for producing potential en- 
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O ' corresponding to the X(5) case, PES with a bump are obtained. The PES corresponding 
to the Pt chain of isotopes suggest a transition from prolate to oblate shapes at ^^^Pt. 



ergy surfaces (PES) for series of isotopes suggested as exhibiting critical point symmetries. 
Relatively flat PES are obtained for nuclei showing the E(5) symmetry, while in nuclei 
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F3 . 

1 Introduction 

X 

' Critical point symmetries |2] corresponding to shape phase transitions in nuclear struc- 
ture are recently receiving considerable attention, since they lead to parameter independent 
(up to overall scale factors) predictions for spectra and B{E2) transition rates which com- 
pare well with experimental data jHl IH 13 EI- The E(5) P and X(5) j2] critical point 
symmetries correspond to special solutions of the Bohr Hamiltonian [7^, in both of which 
an infinite square well potential in the quadrupole (/?) degree of freedom is assumed. In 
E(5), corresponding to the transition from vibrational [U(5)] to 7-unstable [0(6)] nuclei, a 
7-independent potential is assumed |H], while in X(5), related to the transition from vibra- 
tional [U(5)] to axially symmetric prolate [SU(3)] nuclei, a potential of the form «(/?) +^(7) 
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is used. A systematic study of phase transitions in nuclear collective models has been given 
in PCnilll]. 

It is of interest to examine if the assumption of a relatively flat potential in the E(5) and 
X(5) models is justified through the use of a completely different method, as the Relativistic 
Mean Field (RMF) approach ^21 using the NL3 effective force, the parameter set of which 
has been fixed by fitting ground state properties of spherical nuclei [13]. The calculation of 
Potential Energy Surfaces (PES) for a series of isotopes in which a critical nucleus appears 
should result in a relatively fiat PES for this particular nucleus. 

Furthermore, the question of a prolate to oblate shape phase transition, expected to 
take place in the Pt region, has been raised recently f^. RMF theory ,12j is a suitable 
tool for locating the prolate to oblate shape transition in the Pt series of isotopes. 

The Relativistic Hartree Bogoliubov model used in the present work is briefly described 
in Section 2. In this variation of the RMF theory, pairing correlations are treated in a 
self consistent manner, since in constrained calculations the correct treatment of pairing 
is essential for a reliable description of the potential landscapes. RMF theory is then 
used to produce potential energy surfaces for chains of isotopes involving nuclei suggested 
to be good examples of the E(5) and X(5) critical point symmetries in Sections 3 and 4 
respectively, while in Section 5 a study of isotopic chains related to the prolate to oblate 
shape transition is given. The main conclusions arising from the present results and plans 
for further work are finally discussed in Section 6. 



2 The Relativistic Hartree Bogoliubov model 

The model describes the nucleus as a system of Dirac nucleons which interact in a relativis- 
tic covariant manner through the exchange of virtual mesons. The relativistic extension 
of the Hartree Fock Bogoliubov theory is described in detail in Refs. [131 116j . The gener- 
alized single-particle Hamiltonian of HFB theory contains two average potentials: the self 
consistent field f which encloses all the long range ph correlations, and the pairing field 
A which sums up the j»2)-correlations. In the Hartree approximation for the self-consistent 
mean field, the Relativistic Hartree-Bogoliubov equations read 

(hr>-m-\ A \(Uu{v)\_ (Uu{v)\ 

\ -A* -hj, + m + Xj [V,ir)) ""Hv^fcWj' 

where m is the nucleon mass and ho is the single- nucleon Dirac Hamiltonian 

-m ■ V + (3{m + g^a) + g^uj'^ + g^r^pl + e ^^ = Si^i, (2) 
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with (T, tu, and p the meson fields, A the electromagnetic potential and Qa-, gw, and Qp the 
corresponding coupling constants of the mesons to the nucleon. 

In Eq. (PJ) the chemical potential A has to be determined by the particle number 
subsidiary condition in order that the expectation value of the particle number operator 
in the ground state equals the number of nucleons. A is the pairing field. The column 
vectors denote the quasi-particle spinors and Ek are the quasi-particle energies. The self- 
consistent solution of the Dirac-Hartree-Bogoliubov integro-differential eigenvalue equations 
and Klein-Gordon equations for the meson fields determines the nuclear ground state. In 
the present variation of the model these equations are solved by expanding the nucleon 
spinors f/fc(r) and Vfc(r), and the meson fields, in terms of the eigenfunctions of a deformed 
axially symmetric oscillator potential. The calculations have been performed using the 
NL3 Lagrangian parametrization ^Hl , which has been proved very successfull in describing 
various nuclear properties at and away from the line of /3-stability. Finally, for the pairing 
field we employ the pairing part of the Gogny interaction . The basic advantage of this 
choice is the fact that the Gogny force has a finite range which automatically guarantees a 
proper cut-off in momentum space. 

3 The E(5) shape phase transition 

The first nucleus to be identified as exhibiting E(5) behaviour was ^^^Ba 0, while ^'^^Pd 
jl8[IT?l] also seems to provide a very good candidate. Further studies on ^^^Ba [20] reinforced 
this conclusion, while ^°^Pd has been suggested |2l| to be an E(5) candidate and ^'^°Pd was 
found 1^ to correspond to a numerical solution using a potential [23 1211 instead of 
an infinite square well potential in the E(5) framework. A systematic search 123 on 
available data on energy levels and B(E2) transition rates also singled out ^^^Xe as a good 
candidate. ^^°Xe has also been suggested [2^] as a possible candidate, in agreement with 
a recent report [2ZI on measurements of El and Ml strengths of ^^^~^^^Xe, which provide 
evidence for a shape phase transition around A ^ 130. 

Potential energy surfaces for 96-ii4p^^ ^^^~^^^Xe, and ^^^~^^^Ba are shown in Figs. 1-3. 
i°opd and ^^^Pd show PES which are quite fiat. The PES of ^^^Pd and "°Pd are also 
quite fiat. In addition, it is clear that ^^^^^^°Xe do exhibit fiat PES, being on the way 
from vibrational behaviour, seen in ^^^Xe (i?4 = 2.044 [2H|), to 7-unstable behaviour in 
the lighter Xe isotopes [22] • Similarly, ^^^~^^^Ba exhibit rather fiat PES, being on the way 
from vibrational behaviour, seen in ^^^Ba {R4 = 2.280 [2Hj)5 to 7-unstable behaviour in the 
lighter Ba isotopes. 
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In general, we remark that the assumption of a rather flat potential in the E(5) critical 
point symmetry is supported quite well by the present calculations in the nuclei which have 
been suggested as good E(5) examples. 

4 The X(5) shape phase transition 

The first nucleus to be identified as exhibiting X(5) behaviour was ^^^Sm followed by 
i5°Nd pUj. Further work on ^^^Sm [SI El ESI El and ^^^Nd [121 ESI ES] reinforced this 
conclusion. The neibouring N=90 isotones ^^^Gd [SHI ISZ| and ^^^Dy [SZl ISH| were also seen 
to provide good X(5) examples, the latter being of less good quality. A systematic study [Hj 
of available experimental data on energy levels and B(E2) transition rates suggested ^^^Ba 
as a possible good candidate, in addition to the N=90 isotones of Nd, Sm, Gd, and Dy. 
^^^Ba has also been recently suggested [39^ as a possible candidate. In parallel, ^^'^Nd has 
been suggested as corresponding to an analytic solution [10] using a potential instead of 
an infinite square well potential within the X(5) framework. 

Potential energy surfaces for ^^''"^^^Nd, ^""^"^^^Sm, ^^^"^^'^Gd, and iso-iss^y shown 
in Figs. 4-7, while PES for ^^'^"^'^'^Ba have already been exhibited in Fig. 3. We remark 
that for all of the above mentioned nuclei [^^^'^^'^Ba, ^^^Nd, ^^^Sm, ^^^Gd, ^^^^Dy] the PES 
is not fiat, exhibiting a deeper minimum in the prolate (/52 > 0) regime and a shallower 
minimum in the oblate (/52 < 0) region. Relatively fiat PES occur for the N=86 nuclei 
i^^Nd, i^^Sm, i^OGd, and ^^^Dy (to a lesser extend). 

The present results for > are in good agreement with the PES obtained for ^^^~^^'^Sm 
by Relativistic Mean Field (RMF) theory in Ref. |3T], using the NL3 force, as well as 
the NLl, NLSH and TMl forces. In Ref. [Hj the PES for ^^^Sm, as well as for neigh- 
bouring nuclei, presents a single minimum in the region /?2 > 0, in agreement with Fig. 

5 . Furthermore the present results for ^"^^"^^^Sm and ^^^Gd are in good agreement with 
Nilsson-Strutinsky-BCS calculations |I2] . 

From the above we conclude that the present results, in agreement with earlier calcu- 
lations, do not predict flat PES for the N=90 isotones, which are the best experimental 
manifestations of X(5). However, the existence of a bump in the PES corresponding to 
good experimental examples of X(5) might be related to the success of the confined /5-soft 
(CBS) rotor model [ISIHI], employing an infinite square well potential displaced from zero, 
as well as to the relevance of Davidson potentials [IHl HO] of the form + Pq/P"^ (where Po 
is the minimum of the potential) in the description of X(5) properties [13 HH]. It can also 
be related to the significant five- dimensional centrifugal effect jlHl, found recently through 
the use of novel techniques allowing for the exact numerical diagonalization of the Bohr 
Hamiltonian [SHlEIliSS- 
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5 The prolate to oblate transition 



The chain of nuclei i^^Hf, i82,i84,i86^^ i88,i90,i92Os, i94,i96p^^ ^^^^ i98,200Hg has been suggested 
|14j . on the basis of experimental data, as exhibiting a transition from prolate to oblate 
shapes at ^^^Pt. The PES shown in Fig. 8 clearly show such a transition, since the prolate 
minimum is lower than the oblate minimum up to ^^^Os, while in ^^^Pt the opposite picture 
appears. The PES for ^^^Pt and especially for ^^^Pt is flat, evolving in ^^^Hg and ^°°Hg 
towards a vibrational shape. The quadrupole moments of the nuclei belonging to this chain, 
shown in Fig. 9, also exhibit the transition from negative values (corresponding to prolate 
shapes) to positive values (corresponding to oblate shapes). It should be noticed that in 
Fig. 9 (as well as in similar subsequent figures) the experimental values should be compared 
to the theoretical results for the total quadrupole moments, while the theoretical results for 
the separate contibutions of the protons and of the neutrons are shown for completeness. 
It is, however, worth remarking that the experimental values are close to the theoretical 
contributions of the protons. 

It should be noticed, however, that the PES of the i84-202p^ isotopes, shown in Fig. 
10, exhibit a transition from prolate to oblate shapes between ^^^Pt (prolate) and ^^'^Pt 
(oblate). In i84-i92p^ minima appear, the prolate minimum being lower than the 
oblate minimum in i84,i86p|.^ with the opposite situation occuring in i88-i92p^ Beyond 
^^^Pt the PES becomes flat, evolving towards a vibrational shape, reached at ^°^Pt. The 
quadrupole moments of i84-202p^^ shown in Fig. 11, also exhibit a transition from prolate 
to oblate behaviour between ^^^Pt and ^^^Pt. 

The prediction that the nucleus ^^^Pt is critical is supported by several pieces of evidence. 
As seen in Fig. 12, the /5i-bandheads (normalized to the energy of the 2f state) exhibit 
a minimum at ^^^Pt, while the crossover of the (normalized to the energy of the 2f state) 
bandheads of the jSi and 71 bands also occurs at the same nucleus. Furthermore, mapping 
the Pt isotopic chain on the symmetry triangle of the Interacting Boson Model 
shows [54J that ^^"^Pt lies very close to the shape phase coexistence region of IBM 156j. 

A similar transition from prolate to oblate shapes is seen in the i88~200Qg Qj^ain, shown 
in Fig. 13. The prolate minimum is lower than the oblate minimum in ^^^^^^^Os, while 
the opposite situation occurs in ^^^'^^^Os. In ^^^Os the PES is already quite fiat, evolving 
towards a vibrational shape in ^°°0s. The quadrupole moments for ^^'^~^°°0s, shown in Fig. 
14, also suggest a transition from prolate to oblate shapes between ^^^Os and ^^^Os. 
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6 Discussion 



The present calculations of Potential Energy Surfaces (PES) in the Relativistic Mean Field 
theory using the NL3 force lead to the following main conclusions. 

a) The assumption made in the E(5) critical point symmetry that the potential in (3 
can be approximated by an infinite square well potential is justified, since rather fiat PES 
are found for nuclei suggested as good examples of E(5). 

b) The PES found for nuclei suggested as good examples of X(5) exhibit a bump, in 
agreement with earlier calculations |H1 021 • The presence of a bump might be related to 
the success of the confined /3-soft (CBS) rotor model [131111, the relevance of Davidson 
potentials in the description of X(5) propetries jUl HHl, as well as to the existence of a 
significant five-dimensional centrifugal effect found recently |l9j . 

c) The PES obtained for the series of Pt isotopes suggest a prolate to oblate shape 
transition at ^^^Pt. 

It is certainly of interest to calculate PES as functions of (32 and 7, in order to examine to 
which extent the above conclusions are influenced by the 7 degree of freedom. In particular, 
one should examine carefully the PES in which two minima appear, since one of the minima 
could be a saddle point. In addition it is of interest to examine if and to which extent the 
present results will be changed if methods beyond the mean field approximation (see, for 
example, Refs. [HZlEHlEni) are employed. 
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Figure 1: (Color online) Potential energy surfaces (PES) for ^^^Pd, calculated using 
Relativistic Hartree Bogoliubov theory with the NL3 force. 
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Figure 6: (Color online) Same as Fig. 1, but for ^^^Gd. 
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Figure 7: (Color online) Same as Fig. 1, but for ^^^Dy. 
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Figure 9: Quadrupole moments for the chain of nuclei of Fig. 8, calculated using the 
Relativistic Hartree Bogoliubov theory with the NL3 force. Experimental data are taken 
fromRef. |2H1- 
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Figure 10: (Color online) Same as Fig. 1, but for ^'^^Pt. 
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Figure 11: Same as Fig. 9, but for ^°^Pt. Experimental data are taken from Ref. [2H]. 
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Figure 12: (Color online) Experimental energies of bandheads of the jSi and 71 bands, 
normalized to the energy of the 2f state of the ground state band, for several Pt isotopes. 
Experimental data are taken from Ref. pHj . 
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Figure 13: (Color online) Same as Fig. 1, but for ^°°0s. 
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Figure 14: Same as Fig. 9, but for ^^^"^^''Os. Experimental data are taken from Ref. [28] . 
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